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ABSTRACT 


A relatively simple method has been developed to optimize the location, 
temperature, and heat dissipation rate of each cooled shield inside an insula- 
tion layer. The method is based on the minimization of the entropy production 
rate which is proportional to the heat leak across the insulation. The re- 
sults show that the maximum number of shields to be used in most practical 
applications is three. However, cooled shields are useful only at low values 
of the overall, cold wall to hot wall absolute temperature ratio. The 
performance of the insulation system is relatively insensitive to deviations 
from the optimum values of temperature and location of the cooling shields. 

Design curves are presented for rapid estimates of the lc.ations and 
temperatures of cooling shields in various types of insulations, and an equa- 
tion is given for calculating the cooling loads for the shields. 



NOMENCLATURE 


i i i 

A Area of heat flow, 

Cp Specific heat of the boiloff vapor, kJ/kg*K 

0 Functional defined by Eq. (14) 

F Functional defined by Eq. (13) 

hfg Latent heat of vaporization of the boiloff liquid, kj/kg 

k Thermal conducti vi ty , W/m-K; with subscripts, coefficients in Eq. (1) 

L Overall thickness of insulation, m* 

m,n Exponents in conductivity function, Eq. (1) 

P T^/T^, Temperature ratio 

q Heat flow rate, W 

R Tc/T h , overall temperature ratio 

s Dimensionless entropy production rate defined by Eq. (5) 

• 

S Entropy production rate, W/K 

t Thickness between walls with single shield between, m* 

T Absolute temperature, K 

x Distance from cold wall, m* 

x * Distance from cold wall in a multi-shield configuration, m* 

X x/t, dimensionless distance* 

X' x'/L, dimensionless distance* 

y Defined by Eq. (8) 

Subscri pts 
C Cold wall 

H Hot wall 

1 i-th shield 

min Minimum 

opt Optimum 

S Shield 

*For systems with single shield L = t, x = x', X = X*. 



INTRODUCTION 


The search for the ultimate, energy efficient insulation system has led 
in the past few years to a fascinating rediscovery and application of some 
fundamental concepts of thermodynamics: specifically, the second law and the 
use of entropy production rates and availability (or exergy) for design 
optimization purposes. The classical approach has been to minimize the heat 
flow between surfaces at different temperatures. 

The concept of a single vapor-cooled shield in an insulation has been 
treated theoretically as far back as 1959 in Scott's classic textbook on cryo- 
genics Cl] and designs employing them were described not much later [2], 
Paivanas, et al., obtained a patent C3] and later reported on the use of uni- 
formly spaced multiple shields which were cooled by the boil-off from the in- 
sulated dewar [4], Eyssa and Okasha [5] considered only radiative heat ex- 
change between shields and minimized the total refrigeration power required. 
Hi 1 a 1 , et al., [6,7] used a similar minimization of refrigeration power as the 
design basis. Related works were reported by Bejan, et al., [B-ll]. 

Recently, Bejan [12] proposed a new point of view, based on the second 
law of thermodynamics , which considers thermal insulations as dissipators of 
useful mechanical power (i.e. the availability or exergy) or, alternately, as 
generators of i rreversibi 1 ity or entropy. Thus, in this method, optimization 
of an insulation corresponds to minimization of either the entropy production 
rate or the irreversibi lity , or the decrease of availability. Various 
applications of this concept to insulation systems have been documented sub- 
sequently [13,14]. 

Our work grew out of an examination of Cunnington's paper [13] who 
utilized a numerical technique to find optimum temperatures at given locations 
for one and two shields for a thermal conductivity function of the form 
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Although several equations seemed to be incorrectly printed we have 
found two of the design curves to be essentially correct. Thus, our purpose 
was 

1. To develop a simple optimization technique; 

2. To generalize the results to a broader class of insulations; and 

3. To develop simple design methods for cooled shields. 

The essentials of this report were already published [15], 
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ANALYSIS 

We accept the previously developed concept that to optimize an insulation 
system is equivalent to minimizing the entropy production rate. In addition, 
we assume one-dimensional heat flow and that the heat capacity of the boil-off 
gas is adequate to do the cooling for all shields and does not impose a re- 
striction on the optimization. In contrast to Rejan [9,11] who has developed 
a constrained optimization Based on the heat capacity of the boiloff we employ 
the argument that in all practical systems the boil-off is generated by 
cooling of some equipment in addition to the heat leakage across the 
insul ation. 

Parallel heat paths, e.g. supports, have not been considered. However, 
each path can be optimized separately using its own thermal conductivity 
function. Then a design decision has to be made whether the two structures 
should be independently cooled at their respective optimum conditions. 

We examine the general situation of an insulation where equivalent ther- 
mal conductivity, k, can be expressed as a two-term function of the absolute 
temperature 

k = k^" 1 + k 2 T n (1) 

where, typically, the first term represents actual conduction with m * 1 and 
the second term represents radiation with n j 3. In the following, m and n 
can be any value except -1. 

The heat flow across a layer of insulation can be expressed in terms of 
Fourier's law 


q dx = Ak dT 


( 2 ) 
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Substituting k from Eq. (1) and integrating across a layer from one end 
at 1, to the other at 2, yields 


x 2 ‘ 


■m + 1 


(T 


.m+1 


- T* +1 ) + 


2 , T n+l 

FTT (T 2 


- t; +1 )]. 


(3) 


Now consider the insulation with a cooled shield at T$ located at x 
between a hot surface at T H and a cold one at Tq, separated by the insulation 
thickness, t, as shown in Fig. la. The entropy production rate for the 
insulation can be determined from the heat flows and temperatures as follows 



where q s = q H - qc * 

The heat flow terms can be expressed in the form of Eq. (3) and the re 
suiting expression can be non-dimensi onal i zed using the following terms 


s 5 |jL where k H . k at r H> 
H 


p=Ia 

■ v 


R = 


^ 2 (m + 1) 
y ' kj (n + l) T H 


n-m 


and 


( 5 ) 

( 6 ) 

(7) 

( 8 ) 




X 



(9) 

d 
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The resulting equation is 

no 

* , [ ,, ([(PR) 1 " * 1 - (PR) 1 " - 1 ♦ (PR)"'] 

+ f[(PR) n+1 - (PR ) n - 1 + (PR) -1 ]} 

+ j {R m [P m+1 - P m - 1 + P* 1 ] 

+ y R n [P n+1 - P n - 1 + P -1 ]} (10) 


Since R, tne overall temperature ratio, is generally known, s is a func- 
tion of P and X, and its extreme value can be found by differentiating it with 
respect to each variable separately and setting the results equal to zero. 
This procedure yields two equations to be solved simultaneously: 3s/3P = 0 and 
3s/3X = U. Because of the regular form of the expressions, one of the final 
two equations contains only a single unknown as follows: 


R m F(m.P) + y R n F(n.P) 

[R" 1 ' 1 D(m,P) + y R n_1 D(n,P] 2 

(ID 
( 12 ) 


- F(m.PR) 4 y F (n ,PR ] 


[0(m,PR) + y D(n,PR)]‘ 


D(m,PR) + Y D(n,PR; 


where the following functionals were used: 
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F(b,B) = R b+1 - B b - 1 + B' 1 (13) 

D(t>,B) = (B + 1) B b - bB b_1 - B" 2 . (14) 


Thus, to find the optimum temperature and location for a shield, Eq. (11) 
can be solved for P, and then X can be calculated from Eq. (12). The heat to 
be removed by the shield, <1$ = <1 h “ * can ^e ^ 0un<3 » as before, from 

Eq. (3). In dimensionless form the equation becomes 




_ 1 - (PR)™* 1 4 y[l 


; pr ) 1 3 


(PrT 1 - R^ 1 4 Y[(PR) n41 - R nn ] 


(lb) 


For multiple shields t; represents the distance between the two surfaces 
surrounainy the i-th shield on either side, T H ^ and Tq ^ are the temperatures 
of these two surfaces, = x^/t; is the location of the shield relative to 
t<, and x! is tne location of the shield relative to the cold wall as shown in 
Fig. lb. To determine the optimum temperatures and locations for multiple 
shields, first we assumed a temperature for the first shield next to the cold 
wall, then we used Eqs. (11) and (12) to find the temperature and location of 
the second shield. This process was repeated for the rest of the shields and 
the hot wall. Thus, each shield was optimized consecutively with respect to 
the two surfaces on either side. With given values of the overall 
temperatuare ratio, R, and of the number of shields, the process requires 


iterative solution. 
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To put the results into proper perspective, the entropy production rates 
can be compared to the thermodynamical ly minimun rate obtainable through spa- 
tially continuous cooling. According to Bejan [12], this rate is 




= -k/ H 

T c 


(tc) 


1/2 T -1 


dT]‘ 


( 16 ) 


This expression was evaluated analytically for the single-term functions 


of k, i.e. for y = 0, and numerically otherwise. 
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RESULTS AND DISCUSSION 

The first set of curves. Figs 2 through 9, show the relative entropy pro- 
duction rates for various thermal conductivity functions and for up to four 
optimally cooled shields as functions of the overall temperature ratio 
R = Tj./T h . The curves show that the entropy production rate increases with 
decreasing values of the temperature ratio, R, and with increasing values of 
the exponent, m and n. Adding shields, of course, reduces the entropy 
production rate; but for most of the practical temperature range, say 
0.01 < R < 0.4, only three shields contribute to significant decreases and 
adding a fourth shield can be considered unnecessary. No shields are useful 
at high values of R; but this "high" range is strongly dependent on the 
exponent of the temperature. The curves developed with k = k^ T 0 *^ for one 
and two shields were very close to those given by Cunnington [13], converted 
appropriately. 

Study of the results of two-term conductivities reveals that the curves 
fall between those obtained for each of the two terms alone. If y is small 
the first term, T™, dominates; whereas if y is large (>10), the second term, 
T 0 , controls. Thus, general conclusions can be drawn from examining the re- 
sults of the single-term conductivities. 

The second set of curves. Figs. 10 through 31, show the optimum 

temperature ratios, T$/T H , and optimin locations, x'/L, of cooled shields as 
functions of the overall temperature ratio, T^/T^, for various thermal 
conductivity functions and with different number of cooled shields. 

Figures 10 and 11 show the optimum single shield temperature ratios, 

PR = T$/T h , and locations, X = x/L, for five conductivity functions. Both of 

these functions generally decrease with decreasing R. The other figures in 

this set show shield temperatures and locations for systems with up to three 



9 


shields and for both single-term and two-term conductivities. The results are 
strongly non-linear. For example, for k^T^ and R = 0.01, the optimum tempera- 
ture ratios for three shields are about 0.09, 0.3, and 0.6 and the optimum 
locations are about 0.05, 0.2, and 0.5. As is to be expected, our uncon- 
strained optimization yields a somewhat better performance per shield than 
Bejan's [9,11] constrained method. 

The sensitivities of the entropy production rates to deviations from the 
optimum values of PR and X are demonstrated in the last set of curves. 
Figs. 32 through 35, for single shields. The sensitivity increases with the 
value of the exponents, m and n, but the curves are relatively flat near the 
minima. A +20 percent change from optimum, for example, has negligible 
effect. Thus, the system is relatively tolerant of deviations from the 

optimum design conditions. 

Calculations with two different conductivities on the two sides of a 
cooled shield show that using the better insulator on both sides always yields 
the optimum condition. However, if for some reason two types of insulations 
have to be used, then the better insulator should be placed on the warm side 


of the shield. 
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(a) 



Hot Wall 




Figure 1 


(b) 

Schematic of the Nomenclature for (a) Single and 
(b) Multiple Shields 



Curve Set 1: Figures 2 through 9 

The effect of optimally cooled shields on 
the entropy production rate for various thermal conducti vities. 
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Curve Set 2: Figures 10 through 31 

Optimal shield temperatures and locations for various thermal 
conductivity functions with different number of shields. 
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Curve Set 3: Figures 32 through 35 

System sensitivity to deviations from the optimum shield 
temperatures and locations for two overall temperature ratios 

with one cooled shield 
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APPENDIX 


COMPUTER PROGRAMS 
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SEPARS and SHIELD 

These two programs are essentially identical, but SEPARS is written in 
PASCAL whereas SHIELD is in BASIC. 

To allow for consecutive calculations of different systems, the program 
always recycles to the starting point. Consequently, the first input 
requested is either a 1, if a calculation is to be performed, or a 0, if no 
more work is to be done. 

Next the program requests input of the insulation's characteristics, 
specifically, the two exponents of the temperatures in the two-term 
conductivity function, the maximum number of cooled shields ( <1 0) to evaluate, 
the value of y, and the temperature ratio of the first shield to the cold 
wall, P(l) = T S1 /T C . The program calculates and presents the characteristics 
of all optimal systems of cooled shields from one shield to the maximum number 
specified in the input. 

The flow chart and a program sample follows. 




























B « B / xim 

t<m . L 4 tn 4 b 

J « J 4 1 



ASS I CN L 4 C 1 3 AND LO 
J * 2 

■ - - -= 
LO > LO / Z 1 C J 3 

J = J 4 1 



ASS I CN LO. L 4 CI 3 AND TCHt I ] 
J ■ 1 

TCHCM « TCHtn / PC J 3 
J * J 4 1 



L 3 C 1 ] * L 4 CI] 

XPLC 1 : ■ It! 3 / L 3 m 
J « 2 

V ' = 

ASSIGN L 3 CJ] AND XPLCJ) 

J * J * 1 



TSHti: . tchcu • pen 
J > 2 

h = 

tshcj: . tshcj-j] * pijj 

J « J 4 1 



or/g/nal page is 
of POOR QUALITY 
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. i 
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SINGLE SPACE IN OUTPUT 

«) 


2 29 

WN 2TE2.N -: ' ') 






230 
2 31 

232 

233 

234 

E®. 

(• 

SINGLE SPICE IN OUTPUT 

») 


FUNCTION PVR! II. E 

REU ) 

REAL. 



135 

var 






234 

A 


IEU. 



13’ 

ItCIN 

(• 

CONFUTE XI • * E •) 



2 31 

1 sC'lNU! ! , 






137 

PVR sEIPUi 






140 

E® 

(• 

CONFUTE IX»‘E •) 



141 







2 42 







143 







2 44 

FUNCTION DIE. II REAL) REAL. 




245 

BED IN 

(• 

PUNCTIOIUL D •) 




M* 0 = ( E« : . C l * FVB ( I X .£)-£/ (PVR (II ,(ll-t)))-U t (SQI(II) ) 

1«7 I®. (• FUNCTIONAL D •) 

:ii 

:4? 

ISO 


58 


1 •) 



59 


131 

131 

131 

134 

153 

134 

137 

.59 

157 

•;40 

u: 

141 

143 

144 

143 

144 
147 
141 
147 
170 
17! 

in 

’.73 

174 

175 
174 
177 
.79 
177 
liC 
111 
112 
113 
294 
195 
: 84 
197 
100 
11? 
199 
i»l 

172 

173 

174 

175 
174 
177 
17! 

l?t 

20C 

20 : 

:o: 

203 

204 
193 
20 . 
107 
209 
20? 
2:0 
111 
212 
113 
214 
123 
2:4 
217 
2:0 
21 ? 
220 
111 
222 

223 

224 

225 
224 
227 
219 
22 ? 


nngTiox HMi uiD.im, 

BEGIN (• FUNCTIONAL f •> 

r «iwiux:.ii*i.i)i)-PW(n.t)-i.»*(i 9 /xd 
ini. <* functional r • ) 
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FUNCTION SIHP50NITCH1 11*1)1 EAL. 
TYPE 

ARR.ARRAY! 1 101) Or REAL. 

VAX 

C.T UR 

ounr . int . 

H . UAL . 

R.l IKTtCXR , 


IXCIN (• CONFUTE HININUN ENTROPY PRODUCTION RATE USING SINPSQN S NUMERICAL INTEGRATION 5CHEKE •> 

DE2TAT .(! O-TtHRl/100 0, 

TOR l .1 TO 101 DO 
BEGIN 

CILI .TCHR.DELTATMl-l); 

TIL: .?VR( (PVRiCIl! .M.CAKA'NP! fRPMPVRlCtll ,NI ) ,9. 51 (CCt! 

END. 

N *T[l)»Ttl9!:. 

TOR A .2 TO 1 00 DO 
BEGIN 

IF Ks i (R D 2 V 2)«2) THEN 
H sH»4 O'YIRI 
ELSE 

H <H«2 9*TIR) 

END. 

5 2 HP SON »!$QR'DEtTAT/3 9*H) ) / ( 2 0.CANA*l(P! /HP! ) 

END. (« CONFUTE NIN1NUB ENTROPY PRODUCTION RATE USING SINPSON S NUNERICAt INTEGRATION 5CHENE •) 


(• NAIN PROGRAM ROOT •) 


BEGIN 

p r ; h 

Read’.n, 
READiPrC). 
VK: le PFC=; 
BEGIN 


00 


(« THIS BIOCI IS USED TO INPUT THE INSULATION THERHAl CONDUCTIVITY. NVKBER •) 
(• Of SHIELDS AND 1ST. SHIELD I COLD VAIL TEKFERATURE RATIO •) 


IKPITH 

REAOLN. 

READ INN. NS .CANA Pill), 

S1NCLESPACE . 

IT CANA. t 0 THEN 

VRITELNC THERNAL CONDUCTIVITY OF THE INSULATION IS I « Rl*T**\H.X 1) 

ELSE 

iECIN 

VRITELN! ’ THERNAL CONDUCTIVITY Or THE INSULATION IS X . Il«T«*\H 1 1/ ♦ I1*T*?\N 1 1), 

VRITELNC tI2MH.l)3/EIl*(N.II)4THOT"(N-H) . \GAJU 7 2) 

END, 

SINCLESPACE. 

S1NCLESPACE. 


NPC.N.l 0. 

NPl «N* 1 9. 

YOB I .1 TO NS DO 
BEGIN 
IN: .1-1. 

INI .1-2. 

BID .9 90999 !. 
CC .9 1; 

DO >1.9. 

COUNT,. 0, 
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ORIGINAL 

0f poor 


PAG* 12. 

Quality 


230 
131 

233 
133 

234 

233 

234 
137 

231 

233 
240 
24: 

242 

243 

244 
145 
244 
147 
240 
247 
250 
15'. 
252 
35! 

234 
255 
254 
357 
250 
25! 
240 
24! 
242 
2 o 3 

244 

245 
244 
24? 
24! 
24- 
2-0 
27: 

272 

273 

274 

275 
2'4 
277 
274 
279 
286 
21 : 
202 
213 

204 

205 
204 
207 
200 
209 
270 
27: 

272 

273 

274 

275 
274 
177 
191 
279 
300 
30: 

302 

303 

304 


<• this uoei ciicuLms im ituitivilt •> 


up tit 

pi .pmutu, 

Hi.pv«(in],H)»r(R.Pti])»cini»pvitim.i)‘r(ii.Pim, 
v2 tsaKPVMm.OM'oiMOoi.PtnuuKA’PVMtn.dM on*D<N.pmi>, 
V3 —5QI(D(H, PI l.CBKi'Dd, Pi ) ) / (FIH.Pl )»CBHA'F(N,Pl ) ) , 

C .(V2/VIWV3, 

Cl »C*DD, 

IT CKO 0 THEM GOTO 100, 
ir ClsO 0 THEN GOTO 200, 

CC-.f-O !I*CC, 

ir ois'CCMO oooooi then goto 200 . 

DC -DO 

100 ru 2 .rid.cc, 

ir (It 1 3 >0 777777 ) OB (1113 (0 .000001 ) THEN 
BEGIN 

in: »im-o 7*cc. 
cc =o i*cc 

END. 

200 COUNT .count. 1 , 

UN*:'. (Cl >0 0) OH UBSCEXO ooooon. 


u -ipvi;in:,(H-i.0)!'O(H.pm).cuu'Pvi!HU],(N-i Di»D(N,pmii, 

i:n: .u/!D(N.pi)»gbha'D<n,pdi, 

i : • : *i!Ei]/(i o ♦ 1 1 1 1 : > . 

v sin/.. pi). caka *f (N. pi i ) / ( i o-itm, 

si i : .vKPWHdm.Hi'HH.pmJtCUU'PVHdtn.Ni'HN.pnin/im, 

SI 1 3 sSIIJKl O.CAHA'NPi /HP 1 1 /HP ; , 

(• IN THIS BLOCK ViBIEBLES IDE BS51CNE0 FOR DIFFERENT SHIELD CONF ICU1ATIONS ») 


:r then 

SECIN 

L2nj.s(i.o-in-m/i:n. 

IT 111 THEN 
IF 1)3 THEN 
BEGIN 
B «: .0, 

L4U3 =0 0, 

TOP. S3 TC IH1 DO 
BECIN 

e sB/iit::. 
i4t::. siom.B 
END. 

L 4 E 1 3 .L4I11M1 O-XtlJM.O, 

10 .I.0-2E12, 

FOR J .2 TO I HI DO 10 .10/11(3), 

L0 «L0/2in. 

I4tll .14111.10, 

TCHII) sim, 

FOR J *1 TO IHl DO TCKI 1 3 .TEHIIJ/PIJl, 

END 

ELSE 

BEGIN 

1413 ) *1 0.(1 0 -SE 12 >*(l .l-XI 13 >/(IE 21 *IE 33 ). 
TCH 133 .R[ 3 mPtlJ*P[lJ) 

END 

ELSE 

BECIN 

lorn .im.d o-imi/nii, 

TCHE 2 J *1122 /PE i 3 
END. 

L3I13 .14(13, 

IPLI11 .im/Lim, 

FOR J *1 TO 1 DO 
BEGIN 

13(21 .13(2-13/11(21, 

2 Pl (21 >IPl( 2 -l 3 . 1 ( 33 / 13 (21 

END. 

TSHIlI.TCHdl'Pdl, 

FOB 3 .1 TO I DO TSHI31 .TSHI3- 1 1*PI3) 

END 

ELSE 
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305 

104 

317 

301 

307 

310 

311 
313 

313 

314 

313 

314 
317 

310 

3 1 1 

330 

331 

333 

533 

334 

335 

534 
337 
330 
539 
330 
33: 
333 

333 

334 

335 
354 
33’ 
534 

339 

340 

341 
343 

343 

344 
3 4'. 
344 
347 

341 
349 

250 

35; 

35; 

3:3 

354 

335 

354 

257 

351 

32? 

340 

34. 

342 

343 

344 
3*5 
344 
347 
141 
349 

370 
37; 

372 

373 

374 

375 
174 
377 

371 
37? 


ICCIN 

urn >i.o. 

TCHU3 *1(3 3 . 
TSHUJ.«TCHm»Pm. 
14(13 >1 0; 

mm. iitij 
nu>. 

TSHIM3 >1.0, 
Tim.rTCHin, 


or»s;kal vg' 

of poor? tv 


sinclespace. 

WIITELNI' NUMBER OF SHIELDS . ',MI; 

WRITELNi NUMBER OF JT1H1TI0MS > ', COUNT 3), 

SINCLESPACE, 

SINCE ESP4CE. 

WRITELNi ’ HEAT REMOVAL ENTtOPT PRODUCTION OPTimiN OPTIMUM' >. 

WIITELNI ' RATE RATE LOCATION TEMPERATURE'), 

WRITELNI' I, 

S LHC L ESP AC E . 


ir ;<u THrN 

TOR J >3 TO I DO TJLJ] .TSHIJ-ll, 
13(1. 13 >E3!!1, 

IIl.lJ >1 0-1(13. 


(• IN THIS BLOCK DIMENSIONLESS HEAT IEM0VAL AND EKTROP! PRODUCTION RATES ') 
(• ARE CALCULATED TOR EACH 5HIELO •) 


fOR J =1 TO I DO 
I EG * N 

Z:V;(PVH(TSH:;..3 l NP:)-PVR(TSHCJ3.HPl>>*l3tJ»tJ/I(j4ll-<pVliTSH2Jl,ll?l)-PVR<TJEJ3.H?l)>*l3tJ:/»IJ]l/K?l 
.UPWRiTSH:j.l’,NPl)-PVRlTSH;J2.IQ'!))4lirJ*13/irj«i;*(PVB(TSHrJ3,NPl).PVS(TlIJ3,NP!)Ml3IJ;yBIJ3)/HP! 
QU: *(2l*CAH»*Z2W(J 4*CAHA»NPl/HPl ! , 

SISHI.il sOIJJ/TSHi.’:. 

WRITELNI ' SHIELD ' , J J, • ‘ :S,QC J 3 9 5 ,' ' ll.SISHIJ] 9.5,' ' . ? , IPLt J] 9 .5 , ' ' S.TSHIJI 9 5), 

END. 


i* riNALLT. OTHER QUANTITIES Or INTEREST ARE CALCULATED IN THIS BLOCK •) 

SINCLESPACE. 

OK'QT :iil O-PWRlTSHdl.MF’ WCAKA-CAHA*PWI<T$Ht:].N?l)>U3UmitMJ*H?!!)m O.CAHA'NPI/HPl ), ' 
OCOLD s IPWR !?SHI 1 1 .KP: ) -PVR ITCH! I ; ,MPl i *CJUU »PWH (TSHI 1 3 ,N?l l-CAHA'PWRlTCHt 1 1 , NPl 1 )*L3( 1 3 Ml( 1 )*MPl ; , 
QCCLD >QCOLDM| 0*CAMA*KPl/MPi) , 

SCOLD .QCCLD iTCHd). 

STOTALCI =SCOLD-QHOT. 

FOR J >1 TO I DO STOTALUI .STOTAD E 1 3 . S I SHC J 1 . 

shin::: .sikpsonitchu n. 

STCTNINH J .ST0TAL- 1 3 /SHINE I ) , 

SKAIC3 =((1 0-PWR (TCH1 1 1 .MTl )»CAMA-CAMA , PVKTCH(I I ,(?1 ) )'( 1 l/TCHtl )-l 0WMPDM1 0*CAHA'KP1/HP1) , 
SNAIINI I I r SKA 1 2 1 2 / SHINE 1 1 . 


IF IONS THEN PI 1*1 3 >1 0 / P 3 . 


SINCLESPACE, 

WRITELNI 

WIITELNI' 

WRITELNI 

WRITELNI 

WRITELN-" 

WIITELNI 

WRITELNI' 

WIITELNi' 

WRITELNI' 

WRITELNI' 

SINCLESPACE. 

SINCLESPACE. 

SINCLESPACE 


COED WALL I NOT WILL TCKPCIATVIE RATIO 

HEAT OUT AT COLD WALL 

HEAT IN AT HOT WALL 

ENTROP7 PRODUCTION RATE AT COLD WALL 

EKTROP! PRODUCTION RATE AT HOT WALL 

MINIMUM ENTtOPT PRODUCTION RATE 

HAIIMUM ENTtOPT PRODUCTION RATE 

TOTAL ENTtOPT PROD RATI WITH \I 2.' SHIELDS 

MA2IHVM / N1NIHUM ENTtOPT PRODUCTION RATIO > 

TOTAL I MINIMUM ENTROPT PIODUCTION RATIO « 


\TCKtIl.l4 4), 

'.OCOLO 144). 

\OHOT 14 4). 

'.SCOLD 14 4). 

’ . -QHOT 14 4); 

'.SHINE II 14 4). 
'.SHAim 14 4). 

. ' , STOTALUI . M .4 1 , 
' . SMAX INC 11.14 4). 

' ,STOTMIN( 11 14:4), 


END. 

?rcH. 

READLH. 

IEADIPFC) 

END 

END 

HOP 
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TO PERFORM COMPUTATION* ENTER 1. OTHERWISE * ENTER 0. 


ENTER > H N NS GAMA PC 13 ' 

WHERE: M 1ST. POWER IN THE THERMAL CONDUCTIVITY EQUATION 

N 2ND. POWER IN THE THERMAL CONDUCTIVITY EQUATION 

NS NUMBER OP SHIELDS 

GAMA — «0 IP USING ONE TERM THERMAL CONDUCTIVITY EQUATION 
>0 IP USING TWO TERM THERMAL CONDUCTIVITY EQUATION 
PC13 -- 1ST. SHIELD / COLD WALL TEMPERATURE RATIO. ALWAYS > 1 

? 1.0 3.0 1 2.5 13.0 

THERMAL CONDUCTIVITY OP THE INSULATION IS K • K1«T»*1.0 + K2*T*t3.0 
CK2*(H+1 >3/CKl*CN+l )3*TH0Tt*CN-M> = 2.30 


NUMBER OP SHIELDS * 1 

NUMBER OP ITERATIONS = 33 


HEAT REMOVAL ENTROPY PRODUCTION 

RATE RATE 


OPTIMUM 

LOCATION 


OPTIMUM 

TEMPERATURE 


SHIELD 


0.43837 


1.83639 


0.36744 


0.23611 


COLD WALL / HOT WALL TEMPERATURE RATIO 
HEAT OUT AT COLD WALL 
HEAT IN AT HOT WALL 

ENTROPY PRODUCTION RATE AT COLD WALL 
ENTROPY PRODUCTION RATE AT HOT WALL 
MINIMUM ENTROPY PRODUCTION RATE 
MAXIMUM ENTROPY PRODUCTION RATE 
TOTAL ENTROPY PROD. RATE WITH 1 SHIELDS 


0.013741 

0.014350 

0.432719 

0.911631 

-0.432719 

1.000503 

18.236148 

2.315503 


MAXIMUM / MINIMUM ENTROPY PRODUCTION RATIO = 18. 226982 

TOTAL / MINIMUM ENTROPY PRODUCTION RATIO = 2.314340 


? 1 


* 


TO PERFORM COMPUTATION. ENTER 1. OTHERWISE. ENTER 0. 


ENTER '■ M N NS GAMA PC ID 


WHERE : 


M 1ST. POWER IN THE THERMAL CONDUCTIVITY EDUATION 

N 2ND. POWER IN THE THERMAL CONDUCTIVITY EQUATION 

NS NUMBER OP SHIELDS 

GAMA — »0 IP USING ONE TERM THERMAL CONDUCTIVITY EQUATION 
>0 IP USING TWO TERM THERMAL CONDUCTIVITY EQUATION 
PC1D — 1ST. SHIELD / COLD WALL TEMPERATURE RATIO. ALWAYS > 1 


1 .0 090 


0.0 23.0 


THERMAL CONDUCTIVITY OP THE INSULATION IS K * K1*T**1.0 


NUMBER OP SHIELDS « 1 

NUMBER OP ITERATIONS * 23 
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HEAT REMOVAL ENTROPY PRODUCTION OPTIMUM 
RATE RATE LOCATION 


SHIELD 1 0.75446 7.03151 


0.35870 


COLD HALL / HOT WALL TEMPERATURE RATIO 
HEAT OUT AT COLD WALL 
HEAT IN AT HOT WALL 

ENTROPY PRODUCTION RATE AT COLD WALL 
ENTROPY PRODUCTION RATE AT HOT WALL 
MINIMUM ENTROPY PRODUCTION RATE 
MAXIMUM ENTROPY PRODUCTION RATE 
TOTAL ENTROPY PROD. RATE WITH 1 SHIELDS 
MAXIMUM / MINIMUM ENTROPY PRODUCTION RATIO 
TOTAL /MINIMUM ENTROPT PRODUCTION RATIO 


0.004293 

0.016030 

0.770687 

3.734070 

-0.770687 

3.504633 

115.966533 

9.994893 

33.089491 

2.851908 


NUMBER OP SHIELDS = 2 

NUMBER OP ITERATIONS = 36 


HEAT REMOVAL 
RATE 


ENTROPY PRODUCTION 
RATE 


OPTIMUM 

LOCATION 


SHIELD 1 0.05470 

SHIELD 2 0.884J1 


2.71297 

4.70678 


0.17445 

0.4B690 


COLD WALL / HOT WALL TEMPERATURE RATIO 
HEAT OUT AT COLD WALL 

HEAT IN AT HOT WALL = 
ENTROPY PRODUCTION RATE AT COLD WALL * 
ENTROPY PRODUCTION RATE AT HOT WALL = 
MINIMUM ENTROFY PRODUCTION RATE * 
MAXIMUM ENTROPY PRODUCTION RATE » 
TOTAL ENTROPY PROD. RATE WITH 2 SHIELDS « 
MAXIMUM / MINIMUM ENTROPY PRODUCTION RATIO = 
TOTAL / MINIMUM ENTROFY PRODUCTION RATIO « 


0.000806 

0.001162 

0.940073 

1.440716 

-0.940073 

3.921467 

619.477774 

7.920388 

157.970919 

2.019751 


TO PERFORM COMPUTATION • ENTER l. OTHERWISE » ENTER 0. 


T 0 

0.175 CP SECS r 12415B CM USED. 
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OPTIMUM 

TEMPERATURE 


0.10732 


OPTIMUM 

temperature 


0.02016 

0.18786 
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PROGRAM SHIELD 

’ 80910 IEH THIS IS A *IAS!C’ H0CRAH TO CALCULATE OPTIMUM TEHPE1ATURES. 

1 10020 OEM tOCITIOKS , DO COOLING LOADS EOS COOLED SHIELDS IN 1 CRYOGENIC 

3 10(31 HEN INSOLATION SYSTEM WHOSE THEliUl CONDUCTIVin TOILOVS THE (ELATION 

4 10040 IEH l*C 1 *T ' NO ♦ Cl *T* N# 

5 00145 REH HCDIEIED IN LATE NOV 1012 
i 00059 1EM 

7 0014 0 IEH DEFINITION OF STOOLS USED 
0 00070 (EH 

» 00000 (EH COLD-SIDE WALL TEMPERATURE TO 

10 00070 (EH NASH-SIDE WALL TEMPERATURE T7 

-.1 00100 REN SPACING BETWEEN SHIELDS AT M AND 1-1 Ll(I) 

12 00119 REM OVERALL THICKNESS Or INSULATION L 

13 00120 REH LOCAL SPACING RATIO, LI ( I ) / tl ( 1-1 ) 12(1) 

14 00130 REM OVERALL SPACING RATIO, 1/11(1). L4il> 

15 00140 (EM 'DISTANCE flOH COLD WALL ) I l LS(t) 

14 00150 REM I-TK SHIELD TEMPERATURE T(I) 

17 30140 REH I-TN SHIELD POSITION RATIO Id) 

11 00;’0 REM I-TH SHIELD TEMPERATURE RATIO P(I) (ALWAYS )!) 

17 00110 REM I-TH COLD-WARM TEMPERATURE RATIO 1(1) (ALWAYS (l) 

20 00170 REM I-TH OIHENSIONLE5S ENTBOPT PRODUCTION RATE S(I) 

21 03175 REM I-TH DIMENSIONLESS HEAT REMOVAL RATE 0(1) 

22 002 10 REM TOTAL DIMENSIONLESS ENTROPY PROD RATE. SKI) 

23 00220 REM MINIMUM ENTROPY PRODUCTION RATE SOU) 

14 00230 REM ENTROPY PROD RATI WITHOUT SHIELDS SV'I) 

21 00240 REM, ENTROPY PROD RATI RATIOS S3=S2/50 AND S4=S7/S0 

24 30250 REM NUMBER Or SHIELDS M 1= OR (10) 

27 00243 REM 

21 33241 DIM •*. < 10 1 3 , Y( 101 ) 

27 00270 PRINT" * 

33 00200 PRINT "INPUT I IP MORE WORK IS TO HE DONE, 0 IE FINISHED’ 

31 00270 INPUT A 

11 90300 if Ait THEN 01350 

3! 00210 PRINT ’INPUT M0 .NO ,H, GAMMA t ?(!>’ 

34 00320 INPUT M0 . NO ,H. CO , P( 1 ) 

35 00325 DIP END t Y !* (M0*I )*Y ‘MO-MO >' (Y‘ ( 1-MO ) ) - 1 / ( Y* Y ) 

24 3032 0 CEE rNE ( Y !x(NO*l)*Y‘NO-NO/(Y‘ 1 1-NO ! ) -1 / !T*Y) 

37 00231 DEE ENE ■' Y < *Y ‘ (HO* 1 )*Y‘ HO* 1 Y- 1 

3! 30340 DEE ENG ( Y i *Y' (NO* 1 '-Y'NO+l /Y-l 

37 00310 PRINT ’ EXPONENT M3s",NC." E2P0NENT NOx’.NO," CAMMAx’.CO 

40 30315 M1*M0*I 

41 0335! NlxNC.i 

42 33340 rCH I = 1 TC M 

43 30370 II*;-; 

4* 00310 12=1-2 

45 00370 R I )= 008001 

44 00400 Cx 1 
47 00410 3x! 

40 0042! PlxFiD'R'I) 

47 0043C Wlx(R 1 1 ) ‘HO'ENF (P( I M«C6*R( I ) * NO* ENC ( P < I > ) ) 

5! 004 35 Wls'. (R(!) ' (MO- 1 1 'END(P ( I ) )*C0*R( I ) ‘ (NO- 1 >«ENE:P( I ) ) ) ’ 2 ) 

5 1 004 3! WJ = -( END(Pi > *C0 *FNC ( P 2 ))‘2>/(FNE(Fl > «C0 • ENC i P 1 ) i 

52 00437 CxWKW:*W3 

53 00440 ClxG'D 

54 00450 IT CKO THEN 00500 

55 00440 IE C!«0 THEN 00570 
54 00470 C.- 1'C 

57 00410 IE ABSU' 00000 ! THEN 00570 

58 00470 C<-C 

57 20503 Sd xRC-.C 

40 00510 IE R(I. > 0 77*777 THEN 00540 

41 00520 IE 1(1) ( I 000001 THEN 00540 

42 00530 COTS 00423 

43 00540 R(!ixR'!)- 7*C 

44 00550 Cx 1*C 

45 00540 GOTO 00420 

44 005’3 Ux-(S(I!-(M0-l!*EW)!P!:))*C0»l(I) , (N0-l)«n(I(P(!))) 

47 00575 11 ( IltUMENDlPl )*C0'EKE(Pl ) ) 

40 30510 l(llxll(!)/(l*l!(l)) 

47 00570 V»(EME(Pl )*C0 , ENC(El ) )/( I-XIl )) 

70 305*5 S(I).V.(R!!!-MO*rRE(P(I))*CI*l(II-BO*fNC(P(l)))/Id) 

71 00574 S(I)»S(I)MI*C0*N1/H|)/H1 

72 99400 ir IM THEN 00470 

73 004 13 L3(D>1 

74 00420 RO(l)xRU) 
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7i 00430 Tl(i).RI(l)*P( 1) 

74 10440 14(1). 1 
77 00450 15(I)«I<!) 

71 00440 SOTO 11130 

77 00470 iJdua-xii-mmn 

00 09400 If 1)1 THEN 00720 

01 09470 2.4(2). 1(1). <l-I(l)l/I(2> 

01 00709 X0(2 ’.R! 2 ) (P( 1 1 

03 00710 SOTO 99730 

04 00720 ir 1)3 THEN 00740 

03 09730 l4(3i.W.!-K:n*(l-I(2)W(I<2)«I;3)) 

04 09749 R0 < 3 ) >R( 3 ) ( (P( 2 )*P (2 ) ) 

07 00719 OCTO 09730 
00 99749 

07 00779 :o;;;.o 
*0 90709 fOR J.2 TO !2 

71 0J’»C B=B(I;iJ) 

72 3009 9 14 ! ! > « L 4 ( 1 1 • B 

73 091:0 NEXT J 

74 00020 I4!!'»l4t!)t(;.|(;)wi 

75 03030 L0.1-I 1) 
h 99049 FDR ;.2 TO II 
7 7 00* 50 L 0 = 1 0 / Z : t J 
70 30169 NEXT < 

•' 090*9 10.19 Till 

ISC : 9 ! 6 :• 14 ( ! ! .14( ! :• «L0 

:o: ooo*o r:-::.r(;/ 

1 92 3J7C9 TO R . : --l TO l: 

193 037 1 9 S3 1 ! i.ROi 1 ) / P . J ) 

144 9072! NEXT J 
195 037 3: 1 3 i 1 « 1 4 ! 1 / 

1 34 03743 1 5 ' 1 ) * X . *. ; 7 L 3 ( 2 ) 

19* 39753 FOR J = 2 TO 1 
1 90 90*60 L 3 f J ) = L 2 9 J - 1 ) / S 2 ; J ) 

107 01=70 :!(J:.15iJ.l)*i:J:/l3<J: 

.19 03*09 NE1T J 
111 9 37 ?: Tl.’.i.RO- J)»P(1) 

..2 3:309 FOR j.2 TO i 
115 313.0 Tl ( J ! iTl i J - 1 < *? ! j 
: 1 4 3 . 32 3 NEXT 
::3 9193: 

;.i Ole*: print* * 

11* 51359 PRINT* 1=" 1 

1 1 0 3104 2 PRINT’ * 

519*9 T2i 1..R0. ;) 

.23 31000 FOR 0.2 TO 1 

121 31979 T2.J<.Ti(J-:« 

122 01,00 NEXT J 
1 2 3 31 1 23 i3:U!).'3!l) 

1 24 3 1 1 0 4 2 1 1 » 1 : 1 1 • I ; 

125 OHIO FOR J=1 TO I 

VIDE LINE 

1 24 31 : 20 21 * U7l ( J*. .' ‘Hi -Tl (J) ‘Hi >*E3( J*1 1 / !< J«2 >-(Tl ( J> ’H!*T2( *H1 )*13(,J) 'X(J) )/H! 

vide ::.sr 

1 27 3 1 1 2 5 Z2«( <T! !J«1 !*Nl*Tl ( J)’K1 )UJ( J*l) /X( J»l)-(Tl ( J> *H1-TX< J) ’HI )UJ!J)/X( J) >iH! 
. 20 0 : 1 20 QlJ):'.2l.CC*22!MUOO*NHHl) 

127 01130 SI ( J ).Q! J ) (Tl ( J ) 

1 30 01 . 49 PRINT* J » " . J . ' Q.*:0(J),* Sl«*.Sl(J>;* 15.M5(J).‘ T/T7«\TliJ) 

131 91150 NEXT J 

1 32 0 1 152 07. C.Tl(!)'R;*C0-CO»Ti(I)*Hl)*mt)/(X(!. 11*111) 

133 01153 Q7*(27/i USO'Nl/Hl > 

1 37 3 1 154 00.(Tl(l!'Hl-R9(:rHl.CO*Tl(i:'Hl-CO*IO(2rHl)*U(lH(XU)*Hl) 

1 35 9 1 155 QO.QD/C.St'NlfHl) 

1 26 91 154 SO.QO/ROd! 

137 01 149 S2M >. SI-07 

1 30 01 142 REH CUCVUTINC OXTi TO SET SHIN 
137 91143 0= i !-R0( 1)1/109 
149 91144 FOR L.l TO 19! 

mi 91145 cai.aouuD'a-i) 

1 42 01 1 44 m)«(C'.l>'RO«GO*Nl/Hl*CU)'KO>'0 J/CU) 

143 03147 NEXT 2 

144 91179 FOR J.l TO 1 

145 01 1 02 32 ( 1 1 .S2 1 2 )«SK Jl 
144 91209 NEXT J 

147 9<29! IEH O0T41N SHIN OSINS SIHPSON S RULE 
MO 91 202 H.T ( 1 >♦!( 10 1 ) 
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14* }!203 
-. 50 41204 

lit ono: 
in ono! 

153 91207 

154 91101 
1!5 01210 
1 34 01220 
127 01230 
HO 01240 

is* oiijc 


49 

01249 

4! 

9127C 

42 

01280 

43 

0 1 2 1 0 

«4 

012P1 

41 

912*2 

44 

512*3 

• 7 

311*5 

4! 

•3 1 2 * ' 

4 * 

51300 

70 

01*51 

? ’ 

31304 

7 * 

0 1 3 1 3 

73 

•31323 

’( 

51 130 

T r 

31340 

7 1 

01350 

5 •) 

1 EOP 


FOR R=2 TO HO 
If l/2*;NTiI/2) THEN 01207 
H.H.2MU) 

CO TO 91200 
a»H*4’f ll 
NEIT I 

S0lI!ti<D/3 , Hl '2 )/ ( 1»C0*N1/H1) 

53(1 )«S2(!l/S0U) 

SK.'is! (!-10!n‘HUCO-CO , ROi!)‘Nli , U/IO(Il-ll/HlWl!«CI*Nl/Hll 

s4(I)«5»<: WS0C! 

IF I.N THEN 01270 
P 1 1 ♦ l ) * l ■' P l 

PRINT’ ’ 

PRINT' P*“ . PI 1 1 , " Rt" ,R( ! ) , “ I=*.I( I! Il.Ml(i), - S:\3w) 
PRINT’ Us" .121 1 ) 14*’ ,14(1) 

PRINT’ ’ 

PRINT’ COLO miim WALL TEKPERATURE RATIO, TO/TP.'JO.'Il 
PRINT' HEAT OUT AT CO ID NALL:" ,Q0 , ’ HEAT IN AT WARN WALL.VQ* 
PRINT* ENTROPT PRODUCTION RATE AT COLD VALL.'.SO 
PRINT" ENTROPY PRODUCTION RATE AT WARN WALL.*,. -OF) 

PRINT* NINIHim ENTROPY PRODUCTION RATE, S0=\50U> 

PRINT’ ENTROPY PROOUCTIOH RATE FOR’ , I , "SHIELDS , Sl.’.SK!) 

PRINT’ HAIimm ENTROPY PRODUCTION 1ATE S*.*,S*(1) 

PRINT" ENTROPY PRODUCTION SATE RATIOS. S3.S2/S0 AND S4=S*/S0" 

PRINT" 53 = " . S3 ( I ) .* S4=’,S4:l) 

NEXT I 
CCTC 00270 
ItC 


OF POOR QL’AuTr 
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NEWRAF 

This program solves the original, complete, constrained optimization 
equations developed in Ref. [9] without the simplifying assumption suggested 
there which eliminated the dimensionless parameter, hf g /C p T H . Only single- 
term thermal conductivity functions were considered in this analysis. 

This program also recycles to the starting point. Consequently, the 
first input is either a 1, if a calculation is to be performed, or a 0, if no 
more work is to be done. 

Next the program requests input of the insulation's characteristics, 
specifically, the exponent of temperature in the thermal conductivity 
function, the number of cooled shields, the dimensionless parameter hf g /C p T H 
for the boiloff from the insulated container, and R = T^/T^. 

The output specifies the optimal characteristics of the given number of 
shields with the constraint that the cooling capacity is limited to the boil- 
off of the liquid due only to the heat leak through the insulation itself. 

The flow chart and a program sample follows. 
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0 0 
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: PROGRAM NEVRAPHI INPUT/, OUTPUT. ARM). 

3 


OTOGWAi. FAGS U 
OF POOR C^ ALl • ‘ 


4 

; 

i 

7 

I 

9 

<0 


*.3 
-;4 
, < 
• M 

■A 

■ 1 ! 

-.9 

:g 

:: 

32 

23 

24 
23 
2 6 

21 
2? 
3C 
•1 
3 2 

ii 

32 

36 


:•! 
•9 
4 5. 

■9 

42 

42 
44 

43 


I. ............. ...... ......... 

(ftl. It. .It. .i..itii.l.tt...l....l..f. I. ... .••••••<•• .................................... Ml) 

(lltl.KI.il ttllt>l.ti.<lt.l>....lt>llttltl ...... l.t..ltttltltttlltlt| III t It l.tll.lll. ...... It ...till. ••.•<. III.) 

.. I. i. .... II. ...It.... ...... III. ..., 

(••llllllllltlltllltlittlltllllltlt ) MEW1AF ( .t... Mill tllllll Mil •«•••• lit Hill) 

(..... .....I. I. .......I... ......... 1........... ....... ...............I, 

(i. .. i. ...... ii. ........ ..... ...... j c CMT0 J j „ iHODADADI I....... .....) 

(. .......... ............ ....III.... DEPT or MECHANICAL ( INDUSTRIAL ENCRC . 

umiv OF ILLINOIS AT URBANA-CHAHPAICH ..••••>•«•«•••»••••••••••• 

............ .. ii......... .. ....... 1204 W CREEK STREET •»..it.«.. 

. UR i ANA, II 41101 


.... ........ ....... ...... .... ..... JULY 1M} i.i...iii..tit.t.t.iii..itit.i.iii., 

{••■t. ...................... ....... ... .a ..... .. ... ............... ii ... j 

............ ............... ............. ................ ........... 

.itittt. titttti.iiitttt. ...i. ..... it ...... .it.i.taitt.tii. 


(t ii it t..iiii.iiiiiiii.i..iiiit iii.iiitiiiiiiititiiii.il t i 
(* ») 
(• THIS PASCAL PROGRAM WAS DEVELOPED TO OPTIMIZE THE ») 
>• LOCATION. TEKPERATURE AND HEAT DISSIPATION RATE •) 
(• or EACH COOLED SHIELD INSIDE AN INSULATION LATER •> 
(* THE THERHAL CONDUCTIVITY OF THE INSULATION HAS *) 
<• THE GENERAL TORN, •) 
(• •) 
<• K . Rl*(T‘*Ni *) 

(» THE OBJECTIVE HAS BEEN TO SOLVE THE SET Of 2*NS.l •) 
(• NON-LINEAR EQUATIONS OBTAINED BT BEJAN. A. "CIS- •> 
<• CRETE COOLING OF LOW HEAT LEAR SUPPORTS TO 42 I," •> 
(» CRYOGENICS. VOL 15, 1975, PP 290-292 •) 
(« *) 
(• SOLUTION 15 BASED ON THE NEVTON-RAPHSON TECHNIQUE •> 
<• DISCUSSED BT STOECKER. W F , DESIGN 01 THERHAL •) 
■ « ST5TEHS.2ND EDITION, SECTION 6-11 , PP 117-119, *1 
(• NCCRAV-H1LL BOOR CO., NT, 1910 *) 
(i i ) 

(.........ii... .............................. .....ii.....) 

(....l..llli■•lll<ltllltl■.l.l.lll■ll•■.l.l.(.ll■ll.ll.l.) 


46 

4? 


41 

2AE22 

; J r. 
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r • 

TTf • 




c < 


ARRAYS. ARRAYS 

i0: or REAL, 

(• 

ii 


ARRAfP. array:: 

111 or REAL, 

(• 

r • 

r * 


arrayt.array:: 

201 or REAL. 

(» 

W N 
EE 

VAR 
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C “2 


A 

ARRAYS, 

(* 

5! 


B 

. ARRAYS, 

(• 

29 


£ 

ARRAYS. 

(• 

60 


D 

. ARRAYS. 

(* 

6 : 


a 

. ARRAYS, 

(• 

6 2 


Ql 

: ARRAYP, 

(• 

6 3 


s 

ARRAYS. 

(• 

64 


SHAI 

: REAL. 

(* 

.5 


SPY 2 K 

. REAL. 

(• 

66 


SHAUN 

- REAL 

(• 

67 


STOTAL 

. RIAL. 

(• 

6! 


STOTMIN 

REAL. 

(• 

69 


SISH 

. ARRAYS. 

<• 

70 


TSH 

ARRAYS. 

(• 

?’ 


7SHC 

. ARRAYS, 

(• 

72 


WORK 

: ARRATT, 

(« 

72 


I 

. ARRAYP. 

(• 

74 


IPL 

ARRAYS. 

(‘ 

75 


2R 

ARRAYS, 

(• 

76 





77 
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AKP. 

: TEXT, 

(• 

79 


BETA 

1EAL . 

(» 


THE SIZE OF ARRAYS DETERMINES THE NAIIHUH NUMBER 0T SHIELDS •) 
THE S12E Of ARRAYP IS N5.1 »i 

THE SIZE OF ARRATT SHOULD BE TWICE THE NUMBER OF SHIELDS •) 


LOWER-OIACOKAL ELEMENTS OF THE TRIDIACONAL HATH 1 1 •> 

DIAGONAL ELEMENTS OF THE TRIDIACONAL MATRIX •) 

UPPER-DIAGONAL ELEMENTS Or THE TRIDIACONAL MATRIX •) 

RIGHT-HAND SIDE OF THE SET OF EQUATIONS DURING ITERATIONS •) 
t-TH OIMENS IONLESS HEAT REMOVAL RATE •) 

DIMENSIONLESS REAT TRANSFER BETWEEN SHIELDS •) 

DIMENSIONLESS ENTROPY PRODUCTION RATI YOI l-TN LAYER •) 

NAIIHUH DIMENSIONLESS ENTROPY PRODUCTION RATE •) 

HINIMUM OlPfZfiS 10NLESS ENTROPY PRODUCTION RATE •) 

SNA! I SHIN •) 

TOTAL DIMENSIONLESS ENTROPY PRODUCTION RATE •) 

STOTAl I SHIN •) 

I-TN DIMENSIONLESS ENTROPY PRODUCTION RATE •) 

I-TH SHIELD / HOT WALL TEMPERATURE RATIO. ALWAYS <1 «) 

GUESSED I-TH SHIELD / ROT WALL TEMPERATURE RATIO. ALWAYS <1 •> 
DUMMY VARIABLES •) 

SPACING IITVEEN NEICHIORINC SHIELDS / INSULATION THICKNESS •) 
DISTANCE TBOPI COLD WALL / INSULATION THICKNESS •> 

EU) / ni-11) •) 


OUTPUT FILE TO IE USED IF DESIRED •> 
PARAMETER DEFINED IN PROCEDURE INPU7H •> 



ORIGINAL PASS Sri 

OF POOR QUALITY 


19 

ROLD 



REAL. 

(• 

DUMMY VARIABLE USED IN SOLVING THE TEIDIiCQKAl MATRIX <> 

ti 

COUNT 



INTEGER, 

(* 

NUMBER OF ITERATIONS NIEDED TO DETERMINE TSKID'S >) 

12 

DELTATC.DEN 


REAL. 

(• 

DUMMY VARIABLES •) 

12 

OIFF.OIFFRAI 


REAL, 

(• 

DUMMY VARIABLES USED IN CRECXINC CONVERGENCE •) 

14 

DIW.DMAX.OMIN 


REAL. 


DUMMY VARIABLES USED IN SOLVING THE TRIDIACOKAL MATRIX •> 

15 

EPS 



REAL. 

(• 

A SMALL VALUE QSED TO OISERVE IF CONVERGENCE IS OBTAINED •> 

14 

ci.cim.cin 


REAL; 

(• 

DUMMY VARIABLES •) 

17 

COLD 



mi. 

l* 

DUMMY VARIABLE USED IN SOLVING THE Tl I DIAGONAL HATH 11 •> 

II 

t v 



INTEGER, 

(• 

INDICES TOR LOOPS •) 

15 

ITER IN 



INTEGER, 

(• 

INDEX USED TO TERMINATE ITERATIONS •) 

»0 

R 



REAL. 

<• 

POWER Or THE THERMAL CONDUCTIVITY EQUATION ») 


NR! 



REAL. 

(« 

EQUALS M-t •) 

>2 

RPI 



REAL, 

(• 

EQUALS M.l •) 

e_i 

NS 



INTEGER, 

(• 

NUMBER OF SHIELDS •) 

*4 

MSP! 



INTEGER . 

(« 

EQUALS NStl •) 

•5 

?rc . 



INTEGER , 

(• 

PROCRAM FLOW CONTROLLER •> 

44 

TCP. 



REAL, 

<* 

COLD WALL / HOT WALL TEMPERATURE RATIO, AlWATS (1 •> 

? " 

Ti.TIRi' 



REAL. 

(• 

DUMMY VARIABLES •> 

»! 

OCOID 



1EAL . 

(• 

BEAT OUT AT COLD WALL •> 

fl 3 

QHCT 



REAL 

(* 

HEAT IN AT HOT WALL •> 

00 

SCOLD 



REAL. 

(• 

ENTROPY PRODUCTION RATI AT COLD WALL •! 

3 i 
02 
02 

itota: 



REAL. 

(• 

SUN OF I ID'S, SHOULD EQUAL 1 AFTER SUCCESSFUL COMPUTATION • 

04 

o; 

34 

PROCEDURE INFUTH. 






o; 

be;:n 

<• 

INPUT OF DATA HEADING 

•) 


06 

wk :teln . 






.’4 

WR CTELK - 

ENTER 


— -i R NS 

BETA TCH 

•JJ 

WRITELN '). 






* 

wr:teln: ‘ 

VKERI 


R POWER IN THE THERMAL CONDUCTIV1TT EQUATION'). 


WRITE IN'- 



NS NUMBER OF 

SHIELDS'). 

• 1 

WRITECM 



BETA - HFC / (CP»7HCT 

'), 

.i 

WRITE in; • 



HFC — 

HEAT Of VAPORIZATION E J/NC3 ' ), 

-.5 

VR.7FI.fr. 



CP .... 

SPEC 

F1C HEAT AT CONSTANT PRESSURE CJ/IC I)'); 

:4 

WRIT ELNV 



TH07 - 

HOT WALL TENPERATURE OKI'); 

•. *4 

WRITELN ; 



TCH — COLD WALL 

/ HOT WALL TEMPERATURE RATIO, ALWAYS < 1'), 

21 

WRITE LN.' •) 






• 3 
1 

* * 

END 

( t 

INF 

T Of DATA HEADING 

*1 


m * 
22 
23 

PROCEDURE RICH. 






34 

Hu .;4 

( t 

PFC 

* ) 



M C 
* - 

vritecn 






2: 

vriteln • 

TO PERFORM COMPUTATION. ENTER I 

OTHERWISE. ENTER O'), 

17 

vsiteln 






23 

eke. 

(• 

PFC 

i! '/ 




’ 3 

30 


2 * 

33' PROCEDURE SINSLEEPfcCE . 

• : <• sincli space in output •) 

34 WRITELNC ') 

35 INI (• SINCIt SPACE IN OVTPUT •) 

34 

37 

3! 

3' FUNCTION PVRUl.E EEAL) REAL, 


4J 

VAR 


4: 

A 

: RIAL, 

42 

BEGIN 

(• COMPUTE XI«»E 

43 

A »E*LN!!2) . 


44 

PWR .EIPlAi 


45 

END. 

(• COMPUTE IX M E 


44 

47 

41 

49 FUNCTION HA10F2(N0l,K02 READ REAL, 

50 BEGIN <• DETERMINES THE LARGEST OF THE TVO GIVEN NUMBERS *1 

53 IF N01ON02 THEN 

52 if no: 'no: then 

53 RAIOF2 >NO! 

54 ELSE 

55 MAIOF2 =N02 

54 ELSE 

57 IUI0F2..N0: 

5* EKE. (» 

59 


DETERMINES THE LUCEST OF THE TVO GIVEN NUK1EIS •> 



ORIGin.-. 

OF POGrt * 

it 

i: 

42 FUNCTION MIN0F2IND1 .N02.REAL) REAL, 

<3 BEGIN (• DETERMINES THE SMALLEST OF THE TVO GIVEN NUMBERS •) 

44 IF N01ON02 THEN 

45 IF N0DN02 THEN 

44 MIN0F2..N02 

47 ELSE 

il HINOF2 =N0 1 

19 ELSE 

70 HINOF2.NO: 

7: END. (• DETERMINES THE SMALLEST OF THE TVO GIVEN NUMBERS •) 

72- 
73 
• ’4 
75 

71 
77 


71 

79 '• MAIN PROGRAM BOD* •) 

89 

8: REG IK 

92 

12 READLN. 

84 READiFFC ! . 

65 WHILE PFC. 1 DC 

Bt BEGIN 

8’ 

88 (• THIS BLOCK IS USED TO INPUT THE INSULATION THERMAL CONDUCTIvm, HUMBER •) 

r (* OF SHIELD:. HFC/ ICP'THCT 1 AND COLD WALL / HOT WALL TEMPERATURE RATIO •) 

90 

INPUT-: 

92 READ LK 

93 READ. M. NS. BETA. TCHi; 

94 SINGLESPACE. 

95 VRITELM* THERMAL CONDUCTIVITY OF THE INSULATION IS K . f 1 *T*» ' ,M 3 1), 

94 WRITELN HFC / (CP'THOT) = '.BETA. 9 5;, 

9 7 SINGLESFACE 

9! SINGLESPACE, 

99 

0; MP; .“!• : 0, 

o: mm: o. 

93 

34 ■* INITIAL GUESSED VALUES FOR TSHIil'S ARE ENTERED •) 

05 

04 DELTATC = C 0-TCH )/(«=*: 0), 

O' FOR TO NS DO TSHGIJI .J'DELTATGfTCK, 

38 

39 {» VARIABLE USED TO CHECK CONVERGENCE CRITERION IS SET AND THE ITERATIVE PROCEDURE •) 

.3 ('OF NEVTON-RAPHSON METHOD IS STARTED •) 


72 


:: EPS =: 0E - : c , 

.3 COUNT .9 
.4 ITEBIN .0. 

.5 REPEAT 

It COUNT .COUNT-1, 

17 FOR 1.1 TO NS DO 

IB BEGIN 

19 . 01 .TSHGIII, 




73 


ORIGINAL PALL Is 


4C 

4! 

CIPl =1.0 
END, 

OF POOR QUALITY 

42 



43 

(• ELEMENTS OF THE TllOIACOHAt Unit ARE COMPUTED •) 


44 



4' 

11!’ .PWRICIP! . HP 1 l-PVKCI ,!!»•< -H'CI ♦HPl » CTCH- ICTA > 1 : 


44 

am =nn , Pvs(c;,m!)*(H*cnii‘(iETi-TCH*cn«c !»(-2 ombeta-tchu), 


4? 

It 1 1 ■ =81 J UNPi ‘PWRIC! ,1041 MCI* < -H* ( BETA-TCH) -C!*(H»2 81), 


41 

C!1 1 iHPl'PWHillP: ,H!*(BETA-TCH*CIH1 ) . 


4? 

DU I »C!H1* l PWRICIP 1 . N?1 )-PVR(Ci ,H?I ) ♦RETA , H?I 'PWRICI , H ) -TCH*dP 1 *PV1 ( C 1 .14 ) *H? 1 *PVE ( C 1 ,K?1>). 

53 

3[l!:.D!!;.C!»<PVR<Ci.H>»(TCH-BETA>-IETA*m»l*PVR<CI,HHTCH*H?l*PVR(CI.H)-Wl*PVR(CI,HPl>); 

21 

DID .D[’3*CIP:*(BETl'PVB(C!P:.Hi-TCH*PWB(CIH.H)l 


52 

EKE. 


23 

lil) =0 0. 


54 

C CMS ! «C.l, 


it 

(• THE TRICIAGONAL MATRIX SOLVER IS SKOVH !K THIS ILCCE 

»> 

s 7 

(• SEE VESTURE , J. R , 1 RAKDB00I Or MVHER1CAL HATS IX 

• ) 

21 

(• INVERSION AND SOLITION CP LINEAR EQUATIONS, SECTION 

• i 

5? 

<• 2 ?, PP 14-35, JOHN VILET 4 SONS, INC , NY, 1941 

• i 

40 




:r s;;:=3 o then core- toe. 


42 . 

BOL2 ' ;B‘ • • 


43 

col: =o: i : / e : : : . 


«4 

vork::: «c:l: 


42 

VORRIK5-:: =8010. 


64 

OHA2 sAsSISIi!!, 


o 7 

CHIN *aes«b:.3>. 


6! 

ros ; =: to ns do 


6? 

BECIN 



div =b:;]-a:!!«eo:o. 



;r c:v=o o then cctc ioc. 


'l 

dhai «HA:or:oHAi,ASSo:v>;. 


1 • 

CHIN =.4!NO:2(OK!N.AES<SIV! , 


*4 

COLO = <0: IC-AC 1 3*C0LD : /DIV. 



vors::: =cold, 


“c 

BOLD =C:!2/0!V 



work : : -ke : =bgld 


" t 

EC 


*5 

TSflLNS: =50LD 


42 

J -Si 


C 7 

for : =: to ns o: 


52 

BES 1 s 


53 

J 5 v* - . . 


54 

colo sVcrr: j: -vop.R; j«ns: *cclo , 


12' 

TS*!J: =C3LD 


16 

sc, 


17 



41 

(« NEWLY CALCULATED VALUES OF TSHUJ’5 ARE COMPUTED ») 


84 



jr. 

FOR ! =1 TO NS DC TSUI!! .TSHSII J-TSHtlC, 


12 

(• CONVERCENCt IS CHECKED IF THE CRITERION IS SATISFIED, 

THE ITEBATION IS •) 

• 3 

! * TERMINATED, OTHERWISE THE NEWLT CALCULATED TSHIIJ'S ARE 

USED AS NEW •) 

14 

(• GUESSES FOR ANOTHER ROUND OF ITERATION 

*> 

14 

DIFFKAl =1 OE-75, 


97 

FOR I : 1 TO NS DO 


98 

BECIN 


99 

DIFF xABSCTSHL ! 3-T5HCC 13 > , 


30 

DIFFP.AX =NAIOi2iDiFF,DIFFHAI) 


o: 

END , 


32 

IF D! FF9UI ! = EPS THEN 


n 

ITEBiN =1 


34 

ELSE 


05 

FOR ! =1 TO NS DO TSHCII; =7SHtn. 


34 

UNTIL ITEBIN:! , 


37 



38 

(• IN THIS HOCK QUANTITIES USED IN DETERHININE THE SHIELDS' SPACIWS UI CONFUTED * 

33 



20 

TOR : =’. TO NS DO 


\ ; 

BECIN 


12 

?: =tsh::i, 


* * 

iP NSOl THEN 


14 

IF KM THEN 


15 

IF IONS THEN 


74 

TIPI sTSHU-11 


17 

ELSE 


11 

TIHI .TSHIM! 


19 

ELSE 




74 


10 

21 

22 

23 

24 

25 
24 
2? 
21 

29 

30 

31 

32 

33 

34 

35 
34 

ii 

31 

■i 

40 

4; 

42 

43 

44 

45 
44 
4’ 
4! 
4» 


54 


P.' 


I 


TIH1 *TCH 
ELSE 

TIK1.TCH; 


OF POOR QL:AL; 


IRtl 2 * ( IIP 1 • PWB < T I .If >• (TI-TIHl ) )/(PWR(TI ,RP1 )-PWS(TIMl ,NP1 > 1 


CKD. 


DON s 1 0. 

rOR 2 >1 TD NS DO DEN .DEN'UCNS-MM .0; 
N5P 1 «KS«: ; 


ri 


(• FINALLY, SPACINCS IETVEEN SHIELDS 2ND OTHER QUANTITIES OF INTEREST ARE CALCULATED •) 

111! <1 0 /DEN , 

IP ill 3 .2113, 
iTOTAi = : : i : . 

FOR : .2 TO MSP 1 DO 
EEC IN 

XE 1 5 «IE I-l 3*XR5 1- 1 : . 

2F iONSPl THEN IPLIIJ «XPIU-13»X[I2; 
rre-TAL ..ITOTALtim 
END. 

IF <AiS(XTDTAl-i 0) ) 1.0E-5) THEN GOTO 100. 

oit : : .ipwr.tsh- :: . up : ) -pwr ctch . up i ; i m x t i 2 «hp i >, 
aitsspi: o-PVR(75Hrns: . mpj 1 > / (Xcmspi j*kp i > . 

FCS I *2 TO .IS DO QL ET 0 .(PWRCSHII I .HPli-PWRCSHIMI .HPA D/dlllWl ) . 
s:ncles?ace. 

vr::e:n. number or shields > ' . ns 2 > . 

WRITELN: NUMBER Or ITERATIONS = '.COUNT 1); 

SINCLESPACE 

sinclssface, 

WR'TELK ' HEAT REMOVAL ENTROPY PRODUCTION OPTIKUH OPTIMUM';. 

WRIT ELN' RATE RATE LOCATION TEMPERATURE ' ) , 

WRITELNi 

SINCLESPACE 
FCR I *• TO *iS DC 
BEGIN 


35 q: :: sOLiNii-QL::;, 

34 SISHi:: sOIII/TSHii: . 

5? WRITELN; SHIELD M.2.' 5,0112 f 3.' ' li.SISHE’J 9:5,' * :P # fPi£!1 9 5,' •rS.TSHIIi 4.51 

5! en: 

19 SINCLESPACE , 

62 SINCLESPACE . 

4. OHCT =ql;ns?:i, 

s: c c cl. : =ql:::. 

c2 SCOLD sQCOLD.'TCH. 

44 STOTAL sSCOLD-OHCT, 

t3 FOR .< =’ TO NS DC STOTAL =STOTAL»SISH[ J3 ; 

4 4 IF Hi )0 0 THEN 

s' SHIN ssORii: O-PWRITCH, (M/2. 0) ) l/IM/2 .0) ) 

40 ELSE 

1 9 SHIN sSQR i LN( O.'TCK)!, 

?C STOTMIN s STOTAL /SKIN, 

'I SHAT = = ( i 0 - PVR (TCH.HPi ! !•( 1 . 0/TCH-5 01/HP1 ) . 

*2 SHAUN .SHAIJ5MIN, 


* t 

SINCLESPACE . 




74 ' 

WRITELN' ' 

COLD 

WALL / HO? WALL TEMPERATURE RATIO 

= ' ,TCH 14 4), 

75 

WRITELN' ' 

HEAT 

OUT AT COLD WALL 

= '.QCOLD 14 4), 

74 

WRITELN! ' 

HEAT 

IN AT HOT WALL 

• ' , QHOT 14 4), 

?• . 

WRITELN! ' 

ENTROPY PRODUCTION RATE AT COLD WALL 

> * .SCOLD 14 4), 

’0 

WRITELN!' 

ENTROPY PRODUCTION RATE AT HOT WALL 

« 1 , -QHOT 14 4); 

79 

WRITELN! ' 

MINIMUM ENTROPY PRODUCTION RATE 

s • .SHIN 14 4). 

00 

WRITELN: ' 

HA1IKUK ENTROPY PRODUCTION RATE 

c ' , SKAX 14 4); 

■ : 

WRITELN' ' 

TOTAL 

ENTROPY PROD RATE WITH ' ,NS 1/ SHIELDS . ' .STOTAL . 14 :4> . 

02 

WRITELN!' 

MAXIMUM / MINIMUM ENTROPY PRODUCTION RATIO 

> \SHAIIN 14 4). 

03 

WRITELN! ' 

TOTAL 

1 MINIMUM ENTROPY PRODUCTION RATIO 

• ' , STOTMIN 14 4), 

04 

100 SINCLESPACE; 




0: 

ir (DlWrO 0) OR 

(1(13* 

0 0) THEN 


14 

DEC IN 




07 

SINCLESPACE. 




00 

WRITELN! ' 

• — 

> CHECI THE ASSEHUY OF COEFFICIENTS TO IE USED IN TRIDIACONAL MATSII 

09 

WIITELNi ' 

... 

) CHECK THE TtiOIiCQNAl HATS II SOLVER 

90 

END. 





IF (ABS( ITOTAI -1 

0) / 

1 0E-5I THEN 


•2 

BEGIN 




93 

SINCLESPACE, 




94 

WRITELN! ' 

... 

> ITOTAI IS NOT IOUAI TO 1.1 <- 


95 

WRITELN! ' 

... 

) COMPUTATIONS ARE NOT CORRECT <- 

-•) 

94 

END. 




97 

SINCLESPACE. 




91 

99 

SINCLESPACE, 





), 

) 



to FFCH , 

01 BUDIN. 

i: luoirrc' 

oo no 

it a: 

15 It OF. 
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ORIGINAL PAGE 13 
°F POOR QUALITY 


TO PERFORM COMPUTATION, ENTER 1. OTHERWISE, ENTER 0. 


T 1 

ENTER > M NS BETA TCH < 

WHERE: M POWER IN THE THERMAL CONDUCTIVITY EQUATION 

NS NUMBER OF SHIELDS 


BETA — HFC / (CPtTHOT > 

HFG --- HEAT OF VAPORIZATION CJ/KDI 

CP SPECIFIC HEAT AT CONSTANT PRESSURE CJ/KG M 

THOT — HOT WALL TEMPERATURE CKO 
TCH COLD WALL / HOT WALL TEMPERATURE RATIO* ALWAYS 1 

1.0 3 0.0143 0.001 

THERMAL CONDUCTIVITY OF THE INSULATION IS K ■ K1»T*#1.0 

HFG / < CPtTHOT > - 0.01450 


NUMBER OF SHIELDS ■ 3 

NUMBER OF ITERATIONS = 9 


HEAT REMOVAL 
RATE 


ENTROPY PRODUCTION 
RATE 


OPTIMUM 

LOCATION 


OPTIMUM 

TEMPERATURE 


SHIELD 1 0.10438 
SHIELD 2 0.23983 
SHIELD 3 0.47781 


I .36143 
1.12395 
0.89782 


0.09719 

0.28870 

0.58568 


0.06685 

0.23076 

0.53219 


COLD UALL / HOT WALL TEMPERATURE RATIO 
HEAT OUT AT COLD WALL 
HEAT IN AT HOT UALL 

ENTROPY PRODUCTION RATE AT COLD WALL 
ENTROPY PRODUCTION RATE AT HOT UALL 
MINIMUM ENTROPY PRODUCTION RATE 
MAXIMUM ENTROPY PRODUCTION RATE 
TOTAL ENTROPY PROD. RATE WITH 3 SHIELDS 
MAXIMUM / MINIMUM ENTROPY PRODUCTION RATIO 
TOTAL / MINIMUM ENTROPY PRODUCTION RATIO 


0.001000 
0.022985 
0.864998 
22.984544 
-0.B6499B 
3.751018 
499.499501 
25.704743 
133.163725 
6 . 852736 


TO PERFORM COMPUTATION, ENTER 1. OTHERWISE, ENTER 0. 


? 1 


ENTER > M NS BETA TCH -I 

WHERE: M POWER IN THE THERMAL CONDUCTIVITY EQUATION 

NS NUMBER OF SHIELDS 

BETA -- HFG / (CP4TH0T) 

HFG HEAT OF VAPORIZATION CJ/KGI 

CP SPECIFIC HEAT AT CONSTANT PRESSURE CJ/KG KI 

THOT — HOT WALL TEMPERATURE CKI 
TCH COLD WALL / HOT UALL TEMPERATURE RATIO, ALWAYS : 1 

1.0 2 0.0154 0.000806 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


THERMAL CONDUCTIVITY OF THE INSULATION IS K » KitTMl.O 

HFG / (CP6TH0T) » 0.01540 


NUMBER OF SHIELDS - 2 

NUMBER OF ITERATIONS - 8 


HEAT REMOVAL 

ENTROPY PRODUCTION 

OPTIMUM 

OPTIMUM 

RATE 

RATE 

LOCATION 

TEMPERATURE 


SHIELD 1 

0.19732 

1.97595 

0.16252 

0.09986 

SHIELD 2 

0.59037 

1.48999 

0.48495 

0.39623 


COLD HALL / HOT HALL TEMPERATURE RATIO - 
HEAT OUT AT COLD HALL = 

HEAT IN AT HOT HALL = 

ENTROPY PRODUCTION RATE AT COLD HALL 
ENTROPY PRODUCTION RATE AT HOT HALL = 

MINIMUM ENTROPY PRODUCTION RATE « 

MAXIMUM ENTROPY PRODUCTION RATE = 

TOTAL ENTROPY PROD. RATE HITH 2 SHIELDS 
MAXIMUM / MINIMUM ENTROPY PRODUCTION RATIO = 
TOTAL / MINIMUM ENTROPY PRODUCTION RATIO ■ 


0.000806 
0.030677 
0.818366 
38.061092 
-0.818366 
3.776103 
619.846992 
40. 708665 
164.149921 
10.780603 


TO PERFORM COMPUTATION. ENTER 1. OTHERWISE. ENTER 0. 

? 0 

0.072 CP SECS. 11471B CM USED. 

/BYE ; 

3KMUFTC COSTS: 255.028 SRUS AT 6.0059 <= 61.50 
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DESINS 

This program optimizes the characteristics of a single cooled shield with 
different insulations on the two sides. Only one-term thermal conductivity 
functions are considered. 

This program also recycles to the starting point; thus the first input is 
1, if a calculation is to be performed, or 0 if no more work is to be done. 

Next inputs are the characteristics of the two insulations, specifically, 
the exponents of temperature in the thermal conductivity functions on the hot 
and cold sides of the shield, a coefficient ratio ALFA (defined in the 
program), the shield to cold wall temperature ratio, P = T^/Tq, and the hot 
wall temperature, T H . 

The output specifies the optimal characteristics of the cooled shield as 
well as other, related information. 

The flow diagram and a program sample follows. 
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ORIGINAL PAGE 19 80 

1 PIOCRAN DirrCOKD(IWWT/.OffT?UT.SIH); OF PCX)R QUALITY 
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) DCS IKS ( •ll•llllllll•tlllt••llllllllltl•lll 
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1 2 
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•J"' 1 
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{ 1 1 1 1 a 



ill* V bluN STitfT naaaaaiaiiaaaaaaaiaaaiiamifmif 

l i 




ulllIU , !l 41101 #•»»*• * *»»• t* •*»» »**ti at tttt t > 

1 7 




M 11 M 11 1 1 1,11,1 s ' 11,11,111 

1 B 

{ mil 



JULY 1913 •tiiiittii»tMiniiia*taaii«fftiata 


(■ini 




ZO 

<iiiii 




2 ; 

* 1,1 




2 2 

(imi 




2 3 

1 * 1 1 ‘ 1 




2! 



(• 

•) 

24 



(• THIS PASCAL 

PROCRAK WAS DEVELOPED TO OPTIKI2E THE •) 

22 



<• LOCATION, TUTPERATURE AND HEAT DISSIPATION IATE •> 

20 



(• FOR A COOLED SHIELD IN A C1T0CENIC INSULATION •) 

20 



(• SYSTER WHOSE THERMAL CONDUCTIVITY HAS THE fORK. •) 

30 



(• 

• ) 

3. 



(• 

( . I1«(T**H) ON THE ROT SIDE. •) 

32 



(• 

R a H»(T»*N) ON THE COLD SIDE •) 

33 



(* 

• I 

34 



(« THE METHOD 

IS BASED OH THE HINIIII CATION Of THE • > 

35 



(• ENTROPY PRODUCTION RATE WHICH IS PROPORTIONAL TO •) 

34 



(• THE HEAT LEAI ACROSS THE INSULATION •) 

* • 



(• 

• ) 

3 B 





J 7 
<0 





42 

iae:: 

100 . 



42 

ikii: 

230. 



44 

LAB £2. 

300. 



44 

47 

var 




40 





4? 


P 

REAL. 

<• SHIELD I COLD VALL TEMPERATURE RATIO. ALWAYS 1 1 ') 

50 


snai: 

REAL. 

(« HAEINUK ENTROPY PIOOUCTIOH RATE RASED OH II*T»*M •) 

5! 


SKA! 3 

REAL, 

(* NAlimm ENTROPY PHODUCTIOH RATE RASED ON R2*T"N ») 

52 


skin: 

REAL. 

(• MINIMUM ENTROPY PRODUCTION RATE BASED ON R1«T»»H *i 

r * 


SKi.NS 

REAL , 

(* MINIMUM ENTROPY PRODUCTION RATE RASED ON I2»T"N •) 

5 4 


STCTAL 

REAL. 

(» TOTAL DIMENSIONLESS ENTROPY PRODUCTION RATE »i 

« t 

y • 


s:sn 

REAL. 

<» DIMENSIONLESS ENTROPY PRODUCTION RATt AT SHIELD •) 

54 


t:h 

REAL. 

(» COLD WALL / HOT WALL TEMPERATURE RATIO. ALWAYS (1 •> 

r • 


TSH 

REAL. 

(« SHIELD 1 HOT WALL TEMPERATURE RATIO. ALWAYS <1 •> 

50 


I 

REAL. 

!• DISTANCE nOH COLD WALL / THICKNESS RATIO •) 

5* 


31 

: REAL. 

<• t 1 (1 l-I) *) 

40 

4! 





42 


cc 

REAL. 

(• DUMMY VARIABLE ') 

43 


COUNT 

. IMTECE1 . 

(• HV3QER Or ITERATIONS HEEDED TO DETERMINE TCH •) 

44 


00 

REAL. 

(• OUW! VARIABLE •) 

45 


Airi 

REAL. 

(• II2MH»l>milMN«l)3 •) 

44 


C. Cl 

REAL, 

(• DUMMY VARIABLES •) 

47 


11 Ci 

. IHTECE1 . 

<« IROEX TO TERMINATE THE SEARCH FOR TCH •) 

to 


K 

REAL. 

(• POWER Or THE THERMAL CONDUCTIVITY ON HOT SIDE •> 

49 


HP; 

REAL. 

(* IQUALS Nil •) 

70 


N 

REAL. 

<* POWER Or THE THERMAL CONDUCT! vm ON COLD SIDE " 

72 


N?1 

REAL. 

(• EQUALS H«t •) 

?2 


rrc 

nmcn. 

(• PROGRAM riOV CONTROLLER •) 

73 


QCOID 

REAL. 

<* HEAT OUT AT COLD WALL •) 

74 


MOT 

REAL. 

(• HEAT IN AT KIT Will *1 

35 


SCOLD 

REAL. 

(» UfTlOPT PIOOUCTIOH RATE AT COLD WALL •> 

74 


SCR 

niT. 

(• OUTPUT THE TO IE USED tr DESIRED •> 

7; 

71 

79 


TSCT 

REAL. 

(* ROT WALL TEMPERATURE III •) 



ORIGINAL PAGE TS 

IC OF POOR QUALITY 

it 

i: FIOCEDUBE IWUTH; 


1) 

ICCIN 

<< iwirr or data heading «i 

It 

VtITELN; 


15 

VUTELNC 

ENTIR • — » If N UFA t TOT <-—•>; 

l< 

niTUir M, 


17 

Wl ITEEN ( ’ 

vhebe « rows of mi thermal conductivity equation on thi rot sior 

II 

VUTILNl' 

8 POVIt or THE THERMAL CONEUCTI VITT EQUATION ON TKl COLS SISt 

II 

VRITElNC 

alfa -- ti2*(«*i)j/m»(*.n J' >, 

10 

VRITtLN! ' 

P SHIELD / COLO VALL TSMPERATU1E IATIO, ALWAYS ID. 

It 

VRITELNC 

THOT - NOT VALL TXH? EtATUBE CUM. 

12 

VlITELNl ' ') 


13 

INI. 

(• INPUT Or DAT! HIUIHC *) 

II 



»5 



II 



• ? 

PROCEDURE P FCH , 


II 

IECIN 

<• PfCH •) 

l» 

vritun. 


'.30 

WRITE IN * ' 

TO PERFORM COMPUTATION, EJCER 1 OTHERWISE . ENTER 1 

13: 

VRITELN 


:C2 

EKE . 

(• PMH •) 

103 



,31 • 



:35 




PRCCEDUR! S 5 NGLE5PACE . 

107 

BEGIN 

t* SINGLE SPACE IN OUTPUT *) 

.01 

VR.'TILN! ' M 


109 

310 

in:. 

! 1 SINGLE SPACE IN OUTPUT *) 

• 4 4 

• • * 

1:3 

FUNCTION PVR II 

E REAP REAL. 

i:< 

PIE 


• 5 « 

A 

REAL. 

’ll 

BEGIN 

<« COMPUTE !!••£ «i 

;;7 

A sE * IN I II '■ . 


:;l 

PVR :EI?1A: 


1 : ’ 

EKE 

<> COMPUTE II * » E •) 

120 
:2: 
: 2 : 



:: 2 

fvkl'ion eie.ii 

REAL; REAL. 


BESIK 

!' FUNCTIONAL C «i 

.25 

C = v E . i 0 f »PVR : II E 1 -E/ :PWR . II . ( t 0-Ein-U l/saillin 

rift 

• • ■ 
lit 

: 30 

ENE . 

(• FUNCTIONAL 0 

FUNCTION E < £ . 1 1 

REAL i REAL. 

; 1 ; 

BEGIN 

(• FUNCTIONAL F •) 

. 2 2 

F * ( PVR : 1 1 . ( E « 

■1 01 ) i-PWRf XI .El-1 0.11 0/11) 

! 1 1 

.39 
: 3 5 

END . 

(• FUNCTIONAL F * ) 

■ i : 

.31 



:3‘ 



no 



.i: 

; i2 

.13 

1 

MAIN PROGRAM BOOT •) 

BEGIN 


1M 

prc.n. 


MS 

HAUH. 


Mi 

bead-, pr; > . 


M? 

VHiLE PFC= ; 30 


Ml 

BEGIN 


M9 



150 

(• 

THIS RLOCI IS USED TO INPUT THE TWO INSULATION THERMAL CONDUCTIVITIES, ») 

15 1 
• 52 
153 

(» 

SHIELD 1 COLO VALL TEMPERATURE RATIO AND HOT VAIL TEWERATURE •) 

INPUTH. 


151 

REACIN . 


155 

READIM.N.ALFA.P .THOT) , 

131 

SINCLESPACE . 


157 

VRITELNC 

THERMAL CONDUCTIVITY OF THE INSOLATION ON THI »T SIDE IS I . U*tT»»\H 1 

151 

VlITELNl ' 

THERMAL CONDUCTIVITT Or THE INSULATION ON THE COLS SIDE IS R • R2‘lT"',N I 

151 

VRITELNC 

(UMR.DmiMtN.D) « '.ALFA I D; 



no 

HI 

142 

143 

m 

HS 

144 

u? 

HI 

147 

170 

171 

172 

173 

174 

175 
174 
177 
171 
::? 
.00 
101 
103 
103 
H4 
105 
.86 
'.07 
: 88 
107 
:«0 
i ? ■ 

t»i 


WITEIHC HOT VALL TEHPEIATVBE . ' ,THOT » 2/ til'), 

SINCLESPACE, 

SINCIESPACE, 


MF1 «H»1 0. 

NPl .«N*1.0, 

TCH >0 000901. 
CC .0 1. . 

DC >1 I. 

COUNT .0, 


ORIGINAL PAGE IS 
OF POOR QUALITY 


(• THIS (LOCK CALCULATES TCH ITEtATIVEL? *) 


1EFEAT 
TSH *MTCH. 

C «0!N,P) *0<N,P) / F C 00 . P > - PVR ( TCH , (2 0-N) )*D(N,T5H) ‘OCR, TSK) /f (H.TSH) /ALfA, 
Cl »G*DD. 

ir CKO. 3 THEN COTO 100. 
ir cue 0 THEN GOTO 200, 

CC .<-0 1)*CC, 

ir ABS'.CCKO 000001 THEN COTO 200, 

OO i-DD . 

100 TCH =?CH»CC, 

If TCH:0 »•»»?» j OS (TCH(0 OOOOOK THEN 
BEGIN 

TCH sTCH-0 7«CC, 

CC >0 l‘CC, 

if aesi c: ) < o oooooi then ins *1 
END. 

200 COUNT =COUNT» 1 , 

UNTIL Cl. 3 0) OS (ABSICCK0 00000 1 ) OB ( I KD= 1 ) , 


173 

174 1? INC--'. THEN 

175 BEGIN 

174 SINGLESPACE. 

177 SINCIESPACE. 

.70 SINCLESPACE. 

177 VRITELNi -•> OPTimm CRITERION CANNOT IE SATISFIED <— '), 

200 VHITELN ' ..... USE SINCLE INSULATION WITH THE LOWER CONDUCTIVITY <-.-'), 

:o: goto ioo 

202 END 

203 

204 (' OTHER QUANTITIES Or INTEREST ARE CONFUTED IN THIS SECTION «) 

205 

204 22 r-ALFA'PVRiTCH. (N-l 0) )*D(N.P) /D(N.TSH) . 

20 " l = 21/(1 0 * 11 ). 

20! QKOT «tl.0-PWR(TSK,NPl) WtCl 0-I)*NTll, 

2 0 7 SCOLD « ALfA* PVR TCH.NPI )MPVR(P .NPl i-1 0) MPWRlTHOT, (N-N) ) , 2'NP! ) , 


210 

SCOLD .QCOLD/TCH 


* I 

STOTAL »<f (H.TSH 

' C I 0-X UALFA'PVR (TCH, N) *f (N,?)lX)/NPl, 

2:2 

SISH = OHCT-OCOLDi TSH, 

1 1 J 

If H=0 0 THE.! 


2H 

spin: =sq:.;ln 

I S /TCH) > 

2:5 

ELSE 


2 I 6 

SHIN: =SQR 

0-PWR(TCH, IK-'Z 05 ) >/ (H/ J D). 

' ' * 

IF N=0 0 THIS 


: : c 

* • 9 

SN1N2 = 5QR : LN( 

f r : : 

I 0 /TCH) ) 

* •) * 

4 • • 

..d: 

SNIN2 rSORKl 

0-PVR(TCH, (N/2 0)))MN/2 111. 

221 

SNAIl =f iN.TCH) /NPl , 

222 

SfiAII =f IK.TCH. /NPl, 

223 



224 



225 

SINCLESPACE, 


224 

VIITELNC 

mvhbir or iterations 

227 

VHITElNc 

COLD VAIL / HOT WALL TINPERATURI RATIO 

221 

VRITELNC 

SHIELD / HOT VALL TEMPERATURE RATIO 

227 

VRITELNC 

SHIELD LOCATION 

230 

VRITELNC 

HEAT OUT AT SHIELD 

231 

VRiTELNC 

HEAT OUT AT COLD WALL 

232 

VRITELNC 

REIT IN AT HOT VAIL 

233 

VRITELN! ' 

ENTROPY PRODUCTION RATI AT COLD VALL 

234 

VRITELNC 

ENTROPY PRODUCTION RATE AT HOT VALL 

235 

VRITELN! 

ENTROPY PRODUCTION RATI AT SHIELD 

234 

VRITELNC 

HIMIXUN ENTROPY PIODUCTION RATI USED ON I1>T”H 

337 

VRITELNC 

MINIMUM ENTROPY PIODUCTIOI RATI BASED OH I2*T»«R 

230 

VRITELNC 

TOTAL ENTIOPY PIODUCTION 1ATE 

237 

VRITELN! ' 

ENTROPY PROD V/O SHIELD BASED ON ll«T*«H 


* ' , COUNT I), 

« ' .TCH 14 4), 

« ' ,TSH 14.4), 

« *.I 14 il, 

< ' , QHOT -QCOLC 144), 

> ‘ , QCOLC 144), 

. ' ,QHOT 144). 

« ' .SCOLD 14 4 > . 

» ' , -OHOT 14 4), 

> ' ,SISH 14 4). 

> ' , SHIN1 14 4), 

« ‘ .SNIN2 14 4). 

* '.ITOTAl 14 4), 

* ' , SNAIl 14 4). 


82 



240 WIITUMC IKTIOPI MOO V/O SHIELD MSEC ON I2*T»*N • ‘.SIUI2 14 4). 

241 200 S1NCLE5PACE, 

242 ' SINCLESPACt. 

243 SINCiESPACE, 

244 PrCH. 

245 IEADLN. 

244 ItADiPrC) 

247 END 

248 END 
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ORIGINAL FAGE £ 
OF POOR QUALITY 


? 


r 


TO PERFORM COMPUTATION. ENTER 1. OTHERWISE > ENTER 0. 


ENTER > M N ALFA P THOT '! 

WHERE: M POWER OF THE THERMAL CONDUCTIVITY EQUATION ON THE HOT SIDE 

N POWER OF THE THERMAL CONDUCTIVITY EQUATION ON THE COLD SIDE 

ALFA — CK2*(M+1)3/CK1*(N+1)3 

P SHIELD / COLD WALL TEMPERATURE RATIO. ALWAYS : 1 

THOT — HOT WALL TEMPERATURE CIO 

1.0 0.0 20.0 4.5 300.0 


THERMAL CONDUCTIVITY OF THE INSULATION ON THE HOT SIDE IS K 
THERMAL CONDUCTIVITY OF THE INSULATION ON THE COLD SIDE IS K 

CK2*(M+1 >3 
HOT WALL TEMPERATURE 


« Kit (Ttt 1.0). 

= K2t ( TttO . 0 > . 

= 20.00 
= 300.00 ck: 


NUMBER OF ITERATIONS 

COLD WALL / HOT WALL TEMPERATURE RATIO 

SHIELD / HOT WALL TEMPERATURE RATIO 

SHIELD LOCATION 

HEAT OUT AT SHIELD 

HEAT OUT AT COLD WALL 

HEAT IN AT HOT WALL 

ENTROPY PRODUCTION RATE AT COLD WALL 

ENTROPY PRODUCTION RATE AT HOT WALL 

ENTROPY PRODUCTION RATE AT SHIELD 

MINIMUM ENTROPY PRODUCTION RATE BASED ON KltTitM 

MINIMUM ENTROPY PRODUCTION RATE BASED ON K2tTl*N 

TOTAL ENTROPY PRODUCTION RATE 

ENTROPY PROD. W/0 SHIELD BASED ON KltTttM 

ENTROPY PROD. W/0 SHIELD BASED ON K2tTitN 
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0.001666 

0.007497 

0.390755 

0.820144 

0.000497 

0.820641 

0.298568 

-0.820641 

109.396253 

3.680131 

40.925828 

178.307751 

299.619216 

598.241762 


TO PERFORM COMPUTATION, ENTER 1. OTHERWISE, ENTER 0. 

■* 0 


0.044 CP SECS, 


10233B CM USED. 
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